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ANNALS OF MATHEMATT 


Vor. III. June, 1887. 


COPLANAR MOTION OF TWO PLANETS, ONE HAVING A ZERO MASS. 
By Dr. G. W. Hitt, Washington, D. C. 


The supposition that two planets circulate about their central body in the 
same plane enables us to dispense with two differential equations of the second 
order in the general problem of three bodies. The further supposition, that the 
mass of one of them is too insignificant to have any sensible effect on the motion 
of the other, enables us to consider the motion of the latter as known and as 
taking place according to the laws of Kepler. Hence, in this case, the two co-or- 
dinates of the planet of zero mass are the only unknowns; and they are given by 
two differential equations of the second order. These suppositions have, approxi- 
mately, place in several cases in the solar system, but I have more especially in 
view the motion of the satellite Hyperion as disturbed by the action of Titan. 
My object in this paper is simply to point out a method of proceeding, which 
may, I think, be advantageously employed in this case. 

Employing the usual notation x, y, 7, for the rectangular co-ordinates and 
radius vector of the planet whose motion is to be determined, 2’, »’, 7’, for the 
corresponding quantities belonging to the acting planet, 7’ the mass of the latter, 
and .V/ the mass of the central body, the differential equations of motion will be 

d*x 2 dy 22 


at 


a’ dt? oy’ 
where 2, the potential function, has the following expression : — 
Q M ef [ I a'x+ yy] 
oe rs KOT ———- = ere per 
v (2? +) V{(e—#’P + (xy —7'F] ‘” 
The co-ordinates of 1’ satisfy the differential equations 
x’ M+m' _, d*y M-+wm , 
oa + 78 a aU, a - aie >a 7 = ©. 
We can, without any loss of generality, assume that the axis of % is directed to- 
ward the lower apsis of #’. Then the integrals of the last-stated differential 


equations are 
x'=a'(cos:’—e’), y =a’ \/(1 —e”) sine’, 
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, 


where ¢’ is derived from the equation 
nt+ci =e’ —e' sine’, 
, ; AM + mm!’ 
a’, ec’, c’ being constants, and w’ being the equivalent of 4] { 7 | , 
° a’ 

It is desirable to know what the differential equations determining + and 
become when expressed in terms of any other variables. For this end Lagrange’s 
canonical form of the equations serves very conveniently. Let the new variables 
be wand s, and employ the subscript (,) to denote the complete differential co- 
efficient with respect to ¢ of any variable to which it is attached. Then 7 stand- 
ing for } (4? + 3°) expressed in terms of 7, s, 7%, 5, Lagrange’s canonical form 
of the equations is 


adoT oT 228 dot oT 28 


dtcu, cu cu’ ates, cs os 
2A Ox D4 
& ce ce 
As we have =a 4,+35-5,4+ : 
cu cs cl 
, 
5 a) a 
o a 4 Oy . oy 
VY = = U - es Ss — bod 
— wa’ a” 2 
we get 
r ay) 2 oy 2 yr ay oy oy) 
‘yr ee Zi 2 -- -- A ” 
i= 5) 4 | a a) ) Jo i e A. + 3 = | MS) 
, cu 7) cu cs cues | 
44 [ eu . 7 4 ‘| i cxroxr ey ey) 
2 “ a sy" T 1 . = at 
2 os \ 3s : ou ot ou o¢) 
oo a oe Je) 3 dy) 2 
(crer oy 24 - 1[ (3) 4 4 ] 
| A - . “ ay . 
Los or" asctj'' * ot J ov 
It is very plain from the form of Lagrange’s equations that if the variables 
~ 


wu and s were so assumed that one of them, # for instance, should disappear at 
once from the expressions for 7 and 2, we should have an integral of the prob- 


A-y a+ 
. @e¢ i ‘ ce a ' 4 ‘ 
lem. For then mid and, integrating, oe constant. This selection, in 
cd ct, ce t) 


a theoretical sense, is always posssible, and in as many essentially distinct ways 

as there are first integrals of the problem, which, in the present case, are four. 

But, although it is easy, in innumerable ways, to make 2 depend on one variable, f 
it is not so easy to make the six factors of the general expression for 7 depend 
solely on the same variable. And, when we inquire what equations must be 
satisfied for this, we find that they are essentially the same as those which are 
satisfied by the Eulerian multipliers. Hence, nothing is gained by approaching 
the problem from this side. 


I propose to take ~ and s so that 


. / Pe ; Sa 
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where y» denotes a function of ¢ supposed known, but, for the present, left inde- 
terminate. From these equations may be derived 


re = PrP (ae + 8°), ala -b ty = pra. 


Hence the potential function, in terms of ~ and s, becomes 


I VM me’ 
2 = —--- 55 - oy — mor ] ‘ 
pr ' +9) y[u—p yey 


In the general expression for 7 we substitute the values 


a a a A 

re , cy , ct , cy P 
== 9X ~- = oY = —=s— 2 

Os Se Oh? Oe OO Ds Pies 


4 ’ 
F's U(vx 
e . d (| ) u 


ay’ op J , 7 
= = yw) I _ Awy') a(n s. 
cl at 


d(yy') _ d(w" 
. - ae ae dt 


The result is 


T= hy?r? (u,? + 5,7) — a? n'y (1 — e”) (us, — si,) 


“— 2a’ ’ adr’ do . do? Pa - 
+ } [ au” ;— I ) vv . % or! p ‘ | » wo" . (2° 2 s*) : 
! r J dt' dt dt? 


For the sake of brevity we may write, /,, /,, 2, being known functions of ¢, 

T = hh, (um? + 5,7) — h, (us, — sa,) + Fy, (0? + 5°). 
This, substituted in Lagrange’s canonical form of the differential equations, gives 
as the equations of the problem, 


d (. 4 Pe ee 
aad) 72 XK er a we’ 
d ads au dh, — oY 
at { 4 4 al 2h, at re at — hss emt Tig 


Let us now adopt a more general independent variable than the time. Call- 
ing this £, let d¢ = dg, in which # may be regarded as a function either of ¢ or £. 
The second supposition will be the more advantageous. In either case as we 
obtain, on integrating, « and s as functions of £, it will be necessary to have the 
values of £ which correspond to given values of the time, and thus the inverse 


function will have to be considered. Then, in terms of the new independent va- 


riable, 
a h, 4 as ( h, ry (AL) 
- — Hh. ——s5= 
a? 3 2) +g +o 
d (h, *) du dy 2 (02) 
— a a 5 — Ahnu = 
Fa 3 az a fits os * 


We can now consider how pg and 4 should be assumed in order that the dif- 
ferential equations may be most simplified. In the first place it appears impor- 
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tant that the potential function 2 should be freed from the independent variable 
£. This is accomplished by putting »= 1. In the second place it seems we 
cannot readily do better than take the eccentric anomaly ¢’ of the attracting 
planet as the independent variable £. Then 
a? r 
dt=-—— de’, and 0 
a’n 


~ 


Also we have 


h,/@=a"n' (1 —e’ cose’), h,=aln'y (1 —e"), th, = al’n' (1 + &’ cose’), 


7 
= | = ; +. M+ gn! _— an! 
pr an 
For the sake of simplicity let the signification of 2 be changed, and, putting 
m’ 





Wiw="* 
2 oe : eae 
“= ae +s) * Ve FF 3) 
Then our differential equations take the form 
q 
a au - & PR 
aE [ — e' cos ¢’) oo | +2y(1 — 0" —(1+e’cose’)u= Sy? 
7 x » as du 20) 
= [« —e'coseé ] —2y(l adit. —(i+e'cose’)s= — 


It will be noticed that the potential function 2 is, by this assumption of va- 
riables, completely freed from co-ordinates expressing the position of the attract- 
ing planet; and that the two factors 1 — e’ cos ¢’ and 1 + e’ cos é’, very simple 
functions of the independent variable ¢’, are the only evidences of the position of 
this body in the differential equations. And, of the four elements of its orbit, ¢’ 
is the only one we have to deal with. 

We propose now to see whether the introduction of elliptic co-ordinates will 
bring about any simplification in the problem. Supposing 

#225, ej, 
ae ae 
_ : ry -~ =I, and 
aQa+rhA Qa, + A, 
be the equations of a confocal ellipse and hyperbola, a, and a, being constants, 


and 4, and A, the new variables destined to take the place of # ands. By elimi- 
nating +,’ from these equations we obtain 





let 


a— a 


=I; 





— 4; 
(a, + A) (@ + A) 
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Pi (@, + 4) (a, + 4, 
whence ay ee \ (« “e 3 1) (¢ ‘ee | . 


; aQ— a& 





The expression of x, in terms of 4, and /, is obtained from this by simply inter- 
changing @, anda,. Thus 


x= [Stet], 


We now proceed to find what 2 becomes in terms of 4, and /,._ By taking 
the sum of the squares of the last two equations we get 


apt+afoat+atAt+h. 
Thus far a, and a, have been left indeterminate, but we now assume 
a, — a4 =}. 
Then w+ s(x, + 4) 4+ 4/7 
= 2a, +4, + A, + 2y/[(@, + A) (@ + 42)) 
=[V(@ +A) + V+ A))’, 
V (2 + 8°) = (a, + A) + (a2 + 4), 
(#—1P +9 =(4,—3/ 4+ x? 
= 2a, + 4 + 4, — 24 [(a@, + A) (Qa, + A)), 
v[(~—1P + 8} =7 (a +4) —V@ + 4), 
u == 2y [(a, + A) (a, + 4)] + §. 
For the sake of brevity we will now put 
V(@+A)=Pf, V (a, + 4) =2. 
Then it is plain 2 may be written 


a. i—¥ y 


p+ pa 
a Anal. de 
p+a' p—9 


We have now to deal with 7. By taking the logarithms of the values of 
x? and «,*, and then differentiating, we obtain 


— (P+ ah + (Pp — 9. 





‘ dx, dh, ; dhy 

ay ath = a+ hy’ 
‘ ax, adh, dh, 

% @th' ati 
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Whence may be derived 


et 4 nd 4 Fad + Lia tar *.. Jar 
4(a1° + d,’) (a, + A) (a, + Ay : aeliiiied T (a, + Ay ‘ 


2) x 
2 ka + 4)(a, + 4) T (a aa x - jp] te 


On substituting in the factor of dA, d/, the values of +,’ and +? it vanishes, and 


the expression takes the form 
Ay § get Ay As er A, 


(a, +) (FA) a) FH) 





4 (ax, — dz3) = = 


Or, in terms of f and g, we have 
adi? + ds? = alee mes 7 ia +7 a dq. 
r= —-_ 


In like manner we get 


uds — sdu =(p + 9) LV AF. dp — \ (t= =4) dq ]. 


The former expression for 7 was 


, if 1s? : oo 
T=4(1 —e’ cose’ aati ad —;, (i—e” SB + $0 + e’ cos ¢’)(w + s°); 
a as 5 





hence, if we abbreviate by ane 
L af 
\ (4 ae i= = 4, 
T= }(1 — e’cos¢’) [a+ é)F + E +51 
° ds" “) ds” 


» 5 d, ad , 2 
a ] (i —ey(p+9)[a% 12 +4(1 +e’ cose’)(p+ QV. 


7 and 2 are somewhat simplified if we adopt variables » and a, such that 


Pret P~-7ee 
Also, for the sake of brevity, put 


batt} =a 4(a—!] =e 
Then we have 


; ‘dew ade do ade 
T=3(1—e'’ cose’) || 14+ 7 | — 2hk |, = 
x (I co ) [ . fa ' = hk ads! =| 


. if ( . 
vane) (a +h ~ $ (1 + e’ cose’) *, 
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By this transformation 2 is considerably simplified; but, as more than off- 
setting this, 7 is rendered complex. As the expression for « in terms of these 


ee == [ (v vy 
—%\ (y ay i]t: 


it will be perceived that / and & are trigonometrical functions of the angles of the 


variables is 


triangle whuse sides are 1, », and 4, which might have been anticipated from geo- 
metrical considerations. Thus it appears no advantage would result from the 
employment of elliptic co-ordinates. 

Returning, therefore, to the quasi-rectangular co-ordinates # and s, it seems 
some advantage would be gained if we adopt a new system of co-ordinates ~ and 
s, such that the new system is expressed, in terms of the old, as follows : — 


“=U - sy (— I), S=U -$y/(— I). 
We can also adopt the trigonometrical exponential corresponding to the arc ¢’ 
as the independent variable. Calling this £ = e¢* ''~'’, an operator Dis adopted, 


, v ' 
equivalent to £ “,, so that D. f' = if" 


ar j 
In terms of the new variables, 2 has the expression 


[_— »y v 


#@ = (us) + A@—n)e—n)~ Sy (a +- s) 
And the differential equations are . 
D {{1— Fel (£4.27) Du) + 2p (i — ee”) Da + [4 be'(8 + 2°) Ja=- s; 
Dili eC + EDs) 2) Ds t [0 CE Es = — 2" 


Only one of these equations need be actually employed, as either can be obtained 
from the other by changing the sign of } (— 1). We have 


oY I— y y 

— 2 * —_— — 4 . >» | L, 
es yu.ys py (4 —1).9/(s — IF 
oY [— yp Lv 

—-2, = F + v 
ou Ve. 7s » (@ — 1). 4 (8 1) 


For the purpose of integrating these equations, we may adopt the method 
of indeterminate coefficients; and we may employ, as proper to represent the 
values of # and s, the infinite series 


' i j h 
wns .6.,¢0°** *". 
s= 2 .8,,;.2 : . 
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= 


Here 7, 7, and # denote positive or negative integers, zero included ; and the sum- 
mation must be extended so as to include all values for 2, 7, or & from — 
to-+%. Theaandc,c’ are constants and functions of the four quantities e’, », 
a and ¢; a and ¢ being two of the four arbitrary constants introduced by integra- 
tion. The two remaining arbitrary constants serve only to complete the two 
elementary arguments which belong to the attracted planet, and, in this method 
of integration, they can pass unnoticed. 

If we suppose that the orbit of the attracting planet is circular, the differen- 


tial equations reduce to the very simple form . 





A oY 

(D+ 1yfvu=>—2 a 
‘ oY 
(D—1fs=>—25 
ou 


And, in this case, an integral can be found. For multiplying the first by Ds, and 
the second by Du, the sum of the equations, thus multiplied, is an exact deriva- 
tive. Integrating, we get 
Du Ds + us + 22=2C, 

C being the arbitrary constant. 

This integral equation may be combined with the differential equations in 
such a way that one of the terms, regarded as the most difficult of expression in 
a developed form, may be eliminated. For example, if this is taken to be the 


- 


term — 4 of 2, the equations serving to determine the a may be 
| [ (2 — I)(s — 1)] 
taken to be 
I—v 
(s—1)D(D+ 2)u4+14D1a Ds + - (us) u+ 3(u—v)(s—v)+ C=0, 
2 y (us) 2 


sy 


— 35+ #(u—v)(s—vy C=0, 
(us)? + 5 (2 )(s ) + 
in which the constant C is not identical with the former C. One of these equa- 
tions suffices, as the other is a consequence of it. The difference of these equa- 
tions is simpler than either of them, and may be of use. It is 


(w — 1) D(D—2)s+4DuDs ++ 
| 


spe -— S$). 

In attempting to derive periodic series for the co-ordinates of Hyperion, it 
appears to me that it will be easier, in the first instance, to assume that Titan 
describes a circular orbit. And, in the next place, to assume that the perturba- 
tions are periodic functions of the mean elongation of the two bodies. And, as 
it may very easily happen that the terms, depending on the second and higher 
powers of the disturbing force, may quite alter the values of the coefficients, it 


D[(«— 1) Ds —(s — 1) Du — 2 (uw —1)(s — 1)) = aa 
' 
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will be well to employ the method of mechanical quadratures. Starting Hype- 
rion from its line of conjunction with Titan, and at right angles to this line, with 
an assumed velocity, trace out its path until the elongation, between the two 
bodies, amounts to 180°, Then, if Hyperion is again moving at right angles to 
its radius vector, the velocity at the start has been rightly assumed. But if not, 
one makes another trial ; and, by interpolating between the two results, a velocity 
is obtained which will more nearly bring about this condition. And continued 
repetition of these trials will enable us to discover, with all desired approximation, 
the velocity which fulfils this condition. When the path of Hyperion, corres- 
ponding to this velocity, has been traced out, it will be easy, by the well-known 
processes of mechanical quadratures, to assign the periodic series representing 
the co-ordinates of the satellite under the supposed conditions. 

When this is done, corrections to the co-ordinates, proportional to the first 
power of the satellite’s proper eccentricity, can be obtained by the integration 
of a linear differential equation. By comparison of these with observation an 
approximate value of this proper eccentricity will be obtained; a thing to be de- 
sired as we seem to know next to nothing about it at present. Also one will be 
enabled to decide whether the motion of the mean anomaly is more rapid than 
that of the mean longitude, as has been asserted, without sufficient reason as it 
seems to me. 

As illustrating this point, suppose that our moon, instead of having an eccen- 
tricity about 0.055, had one about 0.001. Then the variation would be the pre- 
vailing inequality, and the moon would appear to be in perigee always about 
syzygies, and in apogee about quadratures. In consequence the perigee would 
appear to retrograde with reference to the sun as fast as the moon advances with 
reference to the same body. And yet the relation between the motion of the 
argument, denominated the mean anomaly, and the motion of the mean longi- 
tude, would be nearly the same as it is at present. But the position of the peri- 
saturnium of Hyperion has been concluded from its observed shortest and longest 
vadu vectores. This is allowable only when the inequality, called the equation 


-of the centre, is the overpowering one. 


After the terms, proportional to the first power of the eccentricity, have been 
obtained, those factored by the second, third, etc., powers, can be derived by in- 
tegrating differential equations of the same character. 

In applying the process of mechanical quadratures to the motion of Hype- 
rion, one will meet the difficulty of the uncertain value of the mass of Titan. 
But this cannot be avoided; an assumption must be made, and the results after- 
wards corrected by comparison with observation. 
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ON LOGARITHMIC ERRORS. 
By Pror. H. A. Howe, Denver, Col. 
Mr. R. S. Woodward has kindly sent me the results of his computations of 
the average errors of interpolated values, dependent upon first differences, where 
the interpolating factor has the values 0.1, 0.2, . . . 0.9. I have compared these 


theoretical values with the values given by a discussion of the 1000 examples 
mentioned in No. 6 of Vol. I of the ANNaLs. The comparison is given below: 





eee | Themces | | Actual errors.) Dzterevee 
O.1 0.320 | 0338 | —oo18 
0.2 0.303 0.288 | + 0.015 
0.3 0.304 0.321 — 0.017 
0.4 0.290 0.268 + 0.022 
0.5 0.333 | 0.324 + 0.009 
0.6 0.290 0.27 + 0.014 
0.7 0.304 0.321 — 0.017 
0.8 0.303 | 0.289 + 0.014 
0.9 0.320 0.347. | —0.027 


| 


The agreement of the theory with the observations is quite close, for the er- 
rors are carried to thousandths of a unit of the last place of the logarithm tables 
employed. It is strange that the differences are all negative (except one) when 
the interpolating factor is odd, and all positive when the factor is even. Further- 
more, those actual errors which should be identical theoretically agree very 
closely. 

It appears that the table given on pp. 126-7 of Vol. I does not represent 
well either the probable or the average errors, for the separate values of x, the 
interpolating factor. Hence the approximate theory developed in my articles, 
while it gives fair results, when general averages are taken, is not to be trusted 
for separate values of +. 

It is to be hoped that Mr. Woodward may find time to prosecute further his 
admirable researches in this direction, and apply his theory to a variety of prob- 


lems. 
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ON THE FREE COOLING OF A HOMOGENEOUS SPHERE, OF INITIAL 
UNIFORM ‘TEMPERATURE, IN A MEDIUM WHICH MAINTAINS A 
CONSTANT SURFACE TEMPERATURE. 
sy Mr. R. S. Woopwarp, Washington, D. C. 


1. The problem of the cooling of a homogeneous sphere initially heated to 
a uniform temperature divides itself naturally into two cases, namely: first, that 
in which the dissipation of heat at the surface of the sphere goes on independ- 
ently of the surrounding medium, and second, that in which the dissipation is 
conditioned by the surrounding medium. For brevity we have applied the phrase 
“free cooling” to designate the first case. Both cases have been stated by Fourier 
and subsequent writers, but their investigations are, so far as we know, confined 
almost wholly to the second and more complex case. Moreover, they have, for 
the most part, contemplated the subject from the standpoint of the pure analyst; 
they have been too much occupied with the difficulties of their splendid analysis 
to give much heed to the needs of the computer. 

The most interesting example for application of the theory of a cooling 
sphere is presented by the earth. The hypothesis most commonly entertained 
is, that the whole mass of the earth was at a certain epoch heated to a uniform 
temperature, and is now slowly cooling by conduction without sensibly heating 
surrounding space. With the ultimate object in view of tracing out the conse- 
quences of this hypothesis, we have sought to express the solutions of the case 
of free cooling and that of conditioned cooling in such terms that the computer 
can, without undue labor, assign the temperature at any point within the sphere 
for any value of the time. 

In the following pages the case of the free cooling subject to the restriction 
of a constant coefficient of diffusion is alone considered. How closely the 
conditions of this restricted case accord with those actually presented by the 
earth is a question which requires examination in its proper place. The present 
enquiry, however, is directed solely to the mathematical treatment of the prob- 
lem, and the conditions are tacitly assumed to apply to the earth. 

2. The law of cooling of a homogeneous sphere, initially heated to a uni- 
form temperature and cooling in a medium which maintains a constant surface 
temperature, must evidently be such as to give the same temperature for all 
points equidistant from the centre of the sphere. 

Let #, be the initial uniform excess of the temperature of the sphere over 
that of the surrounding medium. Let be the excess of the temperature of any 
shell of radius ry and thickness dr over the temperature of the surrounding me- 
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dium at any time /, after the initial epoch. Denote the radius of the external 
surface of the sphere by 7”; and let the coefficient of diffusion which is here 
assumed to be constant, be symbolized by a’. Then the function of %,, 7, 7, 7, 4, 
and a? which expresses the law of cooling must be such as to satisfy the follow- 
ing partial differential equation : — 
O(ru) 


a? 
ale 
a* ( 
a 
e 


ae | (1) 


oa” 
If we write the functional relation which is to be discovered in the form 
ru wn f (te, %, 7, t, @), 
we must also have ru=0O for r=”, (2) 
rvu=o for r=0, (3) 
ru=ru, for ¢=0, (4) 
ra=O for t=~a, (5) 
Without going through the steps essential to exclude other functions, we 
observe that conditions (1) and (2) are satisfied by the expression. 


— a2 inz/r.\2Z_- 
Ce O° At 10)" sin ux x] 7, 


C, being any constant, and 7 any integer. Hence, since every term of this form 
will satisfy (1) and (2), we may write 
v2 Z- - , 
ma S Ce Url ro)” sin uz rr. (6) 
. 2 I 
When ¢ = 0, (6) becomes in accordance with (4), 
, 


2 DD 
ru, S C,sinuzr!r, (7) 
7 I 


To fulfill this condition we must have 


} 


. 





; y an x : 
1, j rsin az -dr=C, { sin? nz — dr. (8) 
. ‘, e 7, 
oO 0 
; 


If we put =f Zz 


rt, . ‘ a 7 
(8) becomes me f Hsin xt dt = C, f sin? 74 dt, 
“0 0 


27%, COS NT 


whence C= —- ——., 
7 n 


2r,u, (—1)'*! 
hs . (9) 


| 1) 
2 ie hi 


d\i 
‘ 





OF INITIAL UNIFORM TEMPERATURE, ETC. Ti 


Substituting this value of C, in (6), we find the following equation, which ex- 
presses the law of cooling sought, in its simplest form : — 








n= 
27) My oy ( — 1)" = —at(nz/r,)*t _. 7 
r= — wa t Sin wz 
ria Td 7 
siege ( 10) 
a* (| 7,)*t , r 
~~ € ( sin 7- q 
y, 
__ —a‘* (27 a* : ye 
2r, uy 2! sin 2 
= < ‘ 
a* (37 o}* . F 
+ he = ? sin 37 - 
“ 


The equation (10) satisfies all of the conditions (1) to (5), because each term 
satisfies them. 
When 7 = 0, (10) gives for the temperature at the centre of the sphere, 
27%, O 2%, My Ari) 
“w= |. - = ; 
z 6 Q 7 dr 





Evaluating this expression, we find for the temperature at the centre of the sphere 


mio 48 at (anrir.)2 at(27 2 
(7 a ( ro)*t 1 ¢ 1? (3 GPO ccs ; iF (11) 


ee v4 
w == au,ie—* 


When ¢=0 this gives 7 =, as may be readily seen by recurring to (10). 
When ¢ > 0 it appears that 7 < %. 

3. The series (10) and (11) express all the circumstances of the cooling 
sphere under the assumed conditions. They are rapidly converging series when 
a (z/7,)?¢is not much less than unity, since the negative exponents in the succes- 
sive terms increase as the squares of the natural numbers. But in the most im- 
portant application of this theory, namely, to the earth, the value of a°(z/7,)°7 is 
very small unless the time, /, is very large. Thus, for the earth, if the British 





* foot and the year are the space and time units, we have in round numbers 
a? ax 400,” 
7, = 21 000 000 feet. 
Hence, in order to make a*(z/7,)’¢ unity, we must have 
(21 000 000)" 


f= 5 : 
4007" 
* See section (o) of Sir William Thompson’s paper “On the Secular Cooling of the Earth.””. Thomp- 
son and Tait’s Natural Philosophy, Vol. I. Part Il. Appendix (D). The coefficient refers to the unit of 


bulk of the substance of the earth’s crust 
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or about 100,000,000,000 years. When applied to the case of the earth, there- 
fore, the series (10) and (11) are too slowly converging for the purposes of com- 
putation except when the time is very great. 

To overcome this practical difficulty and obtain an expression which will 
readily give the temperature ~, of any point of the sphere at any time, we pro- 
ceed to transform the series (10.) 

4. In the transformation which follows, and again in the sequel, we shall 
make use of a certain Eulerian integral or Gamma function whose value may 
here be briefly demonstrated, the reader being referred for a fuller investigation 
to the better treatises on the integral calculus. The integral is 


DH 
. , 2 _ i 2 
atx z —(13/a) 
py = f e dx cos fx=!"e ** 
. 20 
Oo 


To prove the equality of the second and third members of this equation, differ- 
entiate with respect to ,7 and then integrate by parts. Thus 


os x 
C7 g atx , , 
= = — { ‘ vax sin 3x 
oa ‘ 
Oo 
1¢ x? ° - x - n 
“"" sin t+ , ff, ~ats - 
an, | ee desis Sees f« "ak COS BX. 
2 ae. 
Oo Oo 
ey o 
‘ ct v 
That is A ee AM 
oF 206° 
A a+ 
or oS tei we 
v 2u°’ 
# 
whence log v7 = — +, + logy, 
4a 
an f 28 fay? 
or = We (4p t) 


To determine the constant 7, we observe that when 3 = o, 


x 





atx? 
v= y= fe dx = ! 





e 24 
(@) 
Therefore ™ 
= —(}3/a)? (4,3/ a)? 24 (° —a?*x? - 
yoo by” and ¢ ‘) = = | e ax COS 3X. 
2 Tv 
; | a 
1e) 





Ae I 
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th 


5. Applying the preceding integral to the 7” term of equation (10), putting 


0 


Pav and 4a =+4,?/(a’P), 
we have 


xn 
—a*(nz/7,,)*t sete? | (r,* | a*t) r r,*x? | (4a27) 
all adil é = , (ve AX COS NZX. 


Hence equation (10) may be written 
LH uD 


2n,-u r,?x? (4a*t)—! ’ 3 
ru = _ ¢ ( dx & (—1)"*'x~'cosnxsinuz—. (12) 
; r 
e 
Oo 


ary (xt) 
n= I 
Now, for brevity, put 
a=y/z and @=7r/1r. 
Then (12) becomes 





oo a) 
27 My * —r,2y? (4a2%x27)-1 ‘ ‘ 
en ae 7 (=?) _— dy X (—1)"t!'x~' cos ny sin xd 
se i] ‘a ; 
’ n= I (13) 
{ 2 = 00 ) 
. + ¥(—1)"t'x-'sinun(y + 8) | 
= — Fo My fe ro*y?(4a?r22) —! dy- n=—I t | 
ax) (77). onal . 
L n=I 
\ : f br it yut P ru, 
te or ev r, — : . 
Z gain, J | im = 
_— iy A 
bees 4z*a°t’ 


and denote the difference of the two trigonometric series in the right-hand mem- 
ber of (13) by Q,. Then, the expression to be evaluated is 
rx 


ran f e~O OL @.- (14) 


oO 


6. We must now determine the value of Q;. For this purpose consider the 
well-known equation, 


1 1— or I 9 
0 ¥( ——* —, = + £ COS ¢ + g COS 29 + Lg cos 3¢ +. 
I—2gcos¢+g° 2 


In this, for ¢, substitute = + y + ¢. Then we have 


+(1—z") 


7 3 22 cos (7 a ¢) + g 





— =gcos(vy+¢)—g*cos2(y+¢)+. 


/ 


Nl = 
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Multiply this by ¢y and integrate between the limits + 4 and — 4. The result is 


+ i 7 


. (1 a go?) ade . 
ac— o 7 = _ st rs a] an ros —_ H 
S " F jumpseaite *t eM OTe Fe O-4 


uo 


n= 


— - —tg*sin2(y +4) +4? sin2(y— 4) 
+ de sin 3(y + #)— 1% sin 3 (y — A) 
= Q, say. (15) 
Comparing the series in (13) and (14) with that in (15), it appears that Q, = Q 
when g = 1. Hence equation (14) becomes 
eS H i 2 / LD 
4 ‘ — p2y2 I— g')a¢g ‘ — p2y2 
ru =4P ie { y—yP ae eS »— P*¥® dy, 
. i J ial 7 “ J1+ 2¢gcos(y+¢)+e¢° a 
t oO t oO (16) 


rt 1. 
The first term in the second member of this equation presents no difficulty. Its 


value is 





But in the second term the element-function of the integral with respect to ¢ is 


zero when g = I, except when cos (y + ¢)=—11. In this case the element- 
function becomes 
I — or? I+ g 
o_ = &—x for ga ti. 
au-gr tz 


The value of this integral with respect to ¢ is therefore the value of the single 


I- - 7 i . 
element —~ ° dg when g= 1. To find its value the most direct process ap- 
7” 


pears to be the following* : — 
Replace cos (y + ¢) by its equivalent 1 — 2sin°}(y+¢). Then, after a 
little reduction, the integral becomes 





Now every element of this integral is zero for g= 1 except that for which 


* The mathematical reader will recognize in the integral to be evaluated one which Poisson and other 
writers have used in proving the law of development of a function in a series of periodic terms. The dif- 
ficulty, if it be such, occurs in another form also in the proofs of that law given by Dirichlet and others. 
I am not certain that the evaluation given in the text will commend itself; but it seems to be more direct 


than, and quite as obvious as, Poisson’s process. 
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y+ ¢ =(2u + 1), # being any integer or zero; and the value of the excep- 
tional element is dependent on 4 only so far as it determines within what limits 
y can make y + ¢ =(2”+ 1)7. We have, therefore, simply to consider what 
value the above expression has when y + ¢ varies between (2” + 1) z — 7 and 
(22 + 1) = + 2,2 being an infinitesimal. The indefinite integral is 


I—g 
arc tan [ “ tand(y+ ‘| ; 
ha 7 





: st ‘ ‘ = # ‘ 
and since, when (1 — g) is infinitesimal, ° tan k(y + ¢) varies from — « 
Té Rx 
~ to + ~ as y + ¢ varies over its range, the proper limits of integration are — « 


and + «. Hence we conclude that 


+ ff 


+ 
j—. g* de , some 3 
4 { ( &) ay = [are tan { © tan$(y + ¢) | ]==. 
J 1+2¢gcos(y+¢)+ 4° ors 4 é, 
eo 


4 
= @ _ 


( It remains to find what values of y + ¢ will give cos (vy + ¢) = — 1, or 
what values of y will make y + ¢ an odd multiple of =. Since ¢ may have any 
value between — 4 and + 4, the range of values of y, which make y + ¢ an odd 
multiple of z, will be from 


(22+ 1)z7—4 to (2n+ 1)7+94, 





n being any integer or zero. That is, y must have values lying between 
z—# and <+4, 


z—# and 37+49, 


w 


z—@ and 57+4, 


vw 


7. Recurring now to equation (16) and attending to the limits of y just de- 
rived we have 


ru= — —Pr (f: Pe" dy + (. ~Pee* dy + 





) , 3 ; 
2ay ¢ 204 f . 
27, ee — 
ru = ru, — —" { f: ds + f. oo ae |. (17) 
] a e e 
vo — a _—_ 
27) t 2a, ¢ 
For brevity, let us put m, = r,/ (2a, 2), 


(18) 





ee ea Pea 
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D 
; i aa 
Then observing that fe" dz=rM,, 
| se 
Oo 


we get the following equivalent forms from (17) :— 


mM, +m 3m, + mM 
27 tly a = 
ru = ru, — —— e~"de+ fe-“ds+..:i, (19) 
! “a a 
mM, — mM Zu —m 
Mm. —m 3, — mm 
2 a 2 
ru = ru, — Hylly ( 1— fer" dz — | oe dz—. ! ) , (20) 
Va va J 
oO Mm, +- m 
My +m jy +m 
2r i - - 
ba haa [ fers de + fern de + 7 a }. (21) 
mM, —m 3, — mM 
My —m ju, — Mm 
2 - y 4 ” sa 
r (uy — u) = Hey { — 2 feo" dz — =f "Ge =. 2 } . (22) 
| te J 4 
Oo ml +- WW 


Equations (19) and (20) give the product of the radius ry and the excess of 
temperature w of any shell of infinitesimal thickness at any time. Equations (21) 
and (22) give the product of the radius and the fall in temperature of the same 
shell during the time ¢ 

8. It will now be well to determine whether equations (19) to (22) satisfy 
the conditions (1) to (5) and also to test their equivalence with (10). 

Differentiating equation (20) or (22) with respect to ¢ we find 








} (ru) 27M, pa (wo—m* 2 (m, — m) 
or = oF 
| Fly cam. — me 9 (3% — Mt) 
=- é 0 a 
|= ct 
i 2 (m, + m) 
—— é . 7 A 
a ct 
a 
r—ri, 
But, from (18), Mm, — m= Ps. 
2a 
a + r 4 
1 ee ee, 
2a 
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o(m,—m) rn—r 
nw tgeeaecal 





Hence = t-*= —1(m,— m)t-', 
cl 4a af 
o(m, + m x _ 
* Pes pS he Lim, + m)t—', 
cl 4c = 
a 
c (vu Vil , ' 
and wf = — — [(m, — m) e~ > ~"* — (Mm, + m) erro te 
cl Ly/2 
+ (30m, — we)e— Ome mh +. . |. (23) 


Differentiating (20) with respect to 7, there results, 











aye © (MM, — Mm) ) 
é . 4 
cy 
Ei . Pe (372, — m2) 
e(vu) 27, j 7! ~~ 
rk ‘ 
P (Mm, + m) 
— ¢ , a 
or 
(+. J 
© (7, — mt) I 
By means of (18) Ss — 
y - Oy 2ay t’ 
© (372, — 12) — 
or 2a, ¢t’ 
o(m, + m I 
(my +m) 
cv 2a, 2 
Therefore 
3 (vm) Villy 5m (ity mm) ~ hy +m)? — ( = 
oy ay (x0) Ll mae sa ¥ inne 


Differentiating the last equation again, we find 


f (m,—m)e~r" 
&* (ru) ru, | — (Mm, + myer 
sen ae = (25) 
nd CPV * | + (3, — m)e~ ve" 
a . . . . . . ° . . 











Dividing (23) by (25), we get ee) =e >) , 


which is equation (1). The equations (19) to (24) therefore fulfill the first re- 
quirement. . 
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If we make + = 7,, equation (1g) becomes 


D 


2rU,, = 
t ru =ru,— - € as 
e 


] 4 


oO 


= Vu, — 7, = O, 


and hence # = 0. Condition (2) is thus satisfied. 
If we make 7 = 0, equation (20) gives 


D 
ae fe 
ru =o — rm {1 —4. fe-" ds] =o, 
jag 


as required by condition (3). 
When ¢ = 0, equation (20) becomes 


| ef 2 
ru = FU, — My.) 1 — —- J e-* az 
| ] see 


@) 
= 7M, , 


or N= H,, 


which is in accordance with condition (4). 


Finally, equations (19) to (22) should satisfy condition (5). But if we make 
t= in (21), the result is 


oO oO 
r(u,—u)= ila [ fe -" ds + fe ak ae . 
] rr e 
oO 


oO 


which is ambiguous. The same ambiguity appertains to equation (22). This 
difficulty will appear also when we seek to determine the temperature at the cen- 
tre of the sphere after an indefinitely great, or infinite, time. But it will be ob- 
served that if we can derive an expression from equations (19) to (22) for the 
temperature at the centre of the sphere, and show this expression to be zero 
when ¢ = ~, then those equations must satisfy (5) 
If we make 7 = 0 in (20), it becomes 


-~- i 


“= 1, — Vu, 
Oo 


But the true value of this expression 1s readily found by differentiation to be 


== Ul, [: — Se oe idee OT ee |. (26) 
Va 


r=0O., 


STWR Bee TNE, 


Now since 7, = 7 2a, ¢, the last equation assumes the indeterminate form when 
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UNIFORM TEMPERATURE, 


¢=«. But from the Eulerian integral of § 4 we have 


> 
2 2 —_ 2 = s I 2 
ee ae = ———_ fe-4 "0 dx COS NX. 
My me 
Oo 
Making ~ successively equal to I, 2, 3,. . . in this expression, and substituting 
the equivalents in (26), we get 
a 
u=u,[I 2 fer me dx (cos x + cos 34 + cospa+...)). 
oe 
Oo 


By reference to the series used in § 6 it is seen that when ¢ = 1, 


ie—s). sr ll Eg") 


: —,=cos*#+cos3%*#+ coss*++...., 
1—2¢gcost+ge" 1+ 2gcosx+¢? 3 4 


Hence we have when ¢ = 1, 


Po) 
. . L(y — 1 _ o* 
u=u,[ 1 — 4 ferdime ax(— : Pod ) 2— aU £)_,| ]-e7) 
rid I— 2g COS 1 -+ I+ 2g cos aT eS 
Oo 


When ¢ is infinite, 7, is infinitesimal, and we may for this value of ¢ restrict the 
upper limit in (27) to an infinitesimal 2, since for all greater values of 2, ¢~‘!/"»” 
vanishes. In this case, also, the second term under the sign of integration in 
(27) may be dropped, by reason of the infinitesimal factor (1 — ¢*). Therefore, 
when ¢ = & and g = 1, (27) becomes 


Y9 | 


7 8 


1 fe7 him" (1 — g*) d (4a) 
“u=u,| 1 —- r+g? J} 7 
(25) 
~ I — 2g cos4 + g”* 


— ae 
SS 





Since 2 may be assumed to be of a lower order than 7,, the value of the factor 
—_ (ji ) 


e~ "©" within the limits of integration is unity. The value of the integral 

then, as shown in $6, is =<. Therefore we conclude that when 7 = Oo and¢‘= ~, 
u=u,[1 —(1/z)z]=0, 

and hence that equations (1g) to (22) satisfy condition (5). 

g. Having established the complete equivalence of equations (1g) to (22) 
and (10), we may now consider some practical features pertaining to their appli- 
cations. In the first place it will be observed that the definite integrals in either 
of equations (19) to (22) form an extremely converging series for all but very 
great values of the time. So converging are these series that, for the case of the 
earth, the first term will suffice when ¢ is as great 100,000,000,000 years. When 
¢ exceeds this time, recourse may be had to the original series (10), which then 
becomes highly converging also, as the following numerical examples will show. 
10. If we confine the time to relatively small values, it will suffice to use the 











=o " 
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first term of (1g) to (20). If, further, we restrict ~ to values which do not differ 
sensibly from 7,, we may write (22) in the form 
Pe r 


2a| t 


2 ceiiiaee 
M,—uU=Mu|i— ae 8S 29) 
¢ / | _f é | (29 


This agrees with an expression which Sir William Thomson has used for the 
fall of temperature, in his paper, Ox the Secular Cooling of the Earth.* 

11. To illustrate the numerical application of the formulz derived, it will be 
well to assume such a value of the time that for the earth both sets of equations, 
i. e. (10) and (11), and (19) to (24), and (26) may be applied. 

As a first example let 

ry=tr, and a(z/7,?t=1. 





The value of ¢ required in the case of the earth to give this last equality is, as we 
have seen in $3, about one hundred thousand million years. With these data 


equation (10) gives 


‘= ma (e _ he Jo ae ‘. 
log No. 

e~' 9.5657055 — 10 + 0.3678794 
e~" 6.09135 —10 4 —0.00004I1 
(e~'—te-*+. . .) 9.§656570 — Io 0.367838 3 

4 0.6020600 

= 0.4971499 
“iu, 9.6705671 — 10 0.468 3463 ; 

hence 4= 0.4683463 %,. 


*Thomsoneand Tait’s Natural Philosophy, Vol. 1. Part II. Appendix (D). 

The law of cooling used by Thomson in this paper is that which applies to the diffusion of tempera- 
ture in an infinite solid on the supposition that at the beginning of the time the temperature had two dif- 
ferent constant values on the two sides of an infinite plane. The formula which expresses this law is 


2) (x2) 
2V £7 
= f en ae, 
po. 
Oo 
in which w,, is the half sum and V the half difference of the initial temperatures ; .v is the distance of any 


point whose temperature is v, from the initial plane, the sign of « being positive for the hotter side and neg- 
ative for the cooler side of the plane; ¢ is the time from the initial epoch and Z is the same as our @?. 











Now when c= 0 the formula gives 7 =v, -+ V, and hence the fall of temperature in any time ¢ is 
] y 
1 r 
2) (x?) 21 (44) 
- > , a > . 
V | Vv + # e 27 12 = viy = . v3 dz } i 
-" 7 J | Vy 
Oo re) 
This agrees with (29) if we make Vu, and +=>7r,—r 
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Again, when a*(z/7,)?¢= 1, 2a, ¢= 27, 7; and when + = }%, equations 
(18) give 

My, = $=, 
m =1i7 

my —m = 7 = 0.785 3082, 

m, + mM = $x = 2.3561945, 

37, — Mm = 3x = 3.9269908, 

3m, + m = fr = 5.4977871. 


From the tables of the integral { « ~* ds given in Oppolzer’s Lehrbuch sur Bahn- 
oO 


bestimmung, Vol. Il. Table X, we find for the first two terms in (21) 


i 


fe ~* ds = 0.2355828.78, 


Ft 7 


¢7 
fe ~** dz = 0,0000000.25. 


i 


The sum of these integrals is 


log 
0.2355829 9.3721438 -— 10, 
4. 0.6020600, 
1 = 0.2485749, 
Uy —U > > 
ae ea 0.§316538 9.7256289 — 10. 
Hence by (21) Uy — U = 0.5316538 %, 
or u = 0.468 3462 m,. 


This agrees with the value derived above from (10) within one unit in the seventh 
place of decimals. It will be observed that the contribution of the second term 
in (21) is less than one unit in the seventh place. The first term of (21) would 
therefore suffice when the time is much greater than 100,000,000,000 years for 
all practical computations relative to the temperature of points within the earth. 

12. As a second illustration, let us compute the temperature at the centre of 
the sphere when 7, of (18) and (26) has the value 

= 


WN, = - =e §. 
2a, 7 





rr? 


This gives f= —,, 
4a 


2 eS 2) ee ee 2 SD. A eee a 4 8 ee a SS eee § ES SS eS See. 
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27 3,009,000,009 years. We also have 


| 


which for the case of the earth gives 
for use in (11), when 7, = I, 


Hence the computation by (11) runs thus :— 


No. log 
Ist term. log « 96377843 — 10 
= 0.9942997 
4 0.6020600 
log e(27)" ——-9,0300240 
ge (dr)? 1.0715785 
+ 0.0848050 e— 4)? 89284215 — 10. 
2d term. log ¢ (47)? 0.03002 
4 0.60206 
log v7 0.63208 
a 4.2863 
— 0.C000517 "abe 5.7137, — 10. 


The third and higher terms are insignificant. The sum of the first two is 
0.0847533. Hence by (11) 
u = 0,1095006 %,. 


According to (26) the computation is as follows :— 


log 
+ 0.3678794 7! 9.5657055 — 10 
+ 0.0001234 e7° 6.09135 —I10 
+ 0.3680028 e~'+e-*+... 9.5658511 — 10 
4 0.6020600 
1 = 0.2485749 
0.8304935 (4/Vz)(e~' + e7"+. . .) 9.9193362 — Io, 


Therefore by (26) “= mM, (1 — 0.8304935) 


= 0.1695065 mM. 


Thus the results derived from (11) and (26) for the temperature at the centre of 
the sphere agree within the limits of precision attainable with 7-place logarithms. 
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SOLUTIONS OF EXERCISES. 





32 
A TRIANGLE /QA is inscribed in the triangle ABC. Determine the ratios in 
which ?, Q, & divide the sides BC, CA, AF in order that dQR, BRP, CPQ may 
be respectively 4, 4, 4 of ABC. [W. M. Thornton.) 
SOLUTION, ; 


Put BP = ax, CO = by, AR=cz. Then the hypothesis gives 


_ 
e we 
—— 
~ <Q 
a 
- 


| 
ee 


—_— 
— 


(I—r) y= 
Elimination gives 3547 — 38% + g = 0. 


Whence we find +x, and by substitution, y and zs. The results are 


ee ASD 
35 4200 45 
in which either value of , (46) may be taken. [H. NV. Draughon.] 
82 


THE major axes of two similar and equal concentric ellipses intersect at 
right angles, and the area common to the two curves is half that of either ellipse. 
Find the eccentricity. [Ormond Stone.] 

SOLUTION. 


Let ACB, A’ CP" be quadrants of the two ellipses, the curves intersecting in D. 

Join CD, and draw DVN parallel to AC, intersecting the circle drawn about C 

with radius CB = CB’ in V/, and CP in J, and join CJZ The circular sector 

CBM is formed from the elliptic sector C/D, by reducing the breadths parallel 

to CA in the ratio 6/a. Hence the area of CBM is the same part of the whole 

circle that CAD is of the whole ellipse. But by hypothesis, C/D is ;'; of the 

whole ellipse. Hence C/M is ;'; of the whole circle, and the angle J/CP is 

therefore equal to 4z. Or since by hypothesis DV = CN, 
MN: DN=b:a; 
6/a=tan tz; 
e? = 1 — tan’ dz = 2)//2 — 2. [ 0. i. Mathiot. | 
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GIVEN one vertex of a rectangle and the ratio of its sides, construct the rect- 
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angle so that the extremities of the diagonal opposite to the given vertex shall 
lie 
1. On two given parallel straight lines, 
2. On two given intersecting straight lines, 
3. One on a given line and the other on a given circle, 
4. On two given circles. [R. D. Bohannan.] 
SOLUTION. 


Let C be the given vertex; AZ the required diagonal; and a, 3 the lines on 
which A, # are to lie. When A describes a, A’, the extremity of CA’ = 2. CA, 
will describe a figure a’, similar and similarly placed to a and x times as large. 
B, the extremity of C4 = CA’, will describe an equal figure 3’, located by turn- 
ing 4’ through a right angle about C. The points of intersection of 3’ with 3 
furnish solutions of the problem. The application of the general method to case 
(4) of the problem is easy. In (4) 3’ is a circle cutting / in two points, each of 
which gives a solution. [ WW. MZ. Thornton.) 

{Mr. Henry Heaton sends substantially the same solution. Mr. E, L. Stab- 
ler suggests the generalized form applicable to other curves. ] 


108 

ConsTRuCT an equilateral triangle of given size inscribed in a given equilat- 

eral hyperbola. [R. H. Graves.] 
SOLUTION. 

Any rectangular hyperbola, passing through the vertices of a triangle, passes 
through its orthocentre. Also the nine-point circle is the locus of the centres of 
rectangular hyperbolas passing through the vertices of a triangle. Therefore the 
centre of the required triangle lies on the given hyperbola, and its inscribed circle 
passes through the centre of the curve. Hence the following construction : — 

Construct a diameter equal to two-thirds of the altitude of the triangle, with 
one end of the diameter as acentre and a radius equal to this diameter describe a 
circle; three of the points of intersection of hyperbola and circle are the vertices 
of the required triangle. (Omit point of intersection at end of diameter, for oth- 
erwise the diameter would bisect 2 chord perpendicularly, which is, 27 general, 
impossible.) If two-thirds of the altitude is less than the transverse axis, the 
construction is obviously impossible. [R. H. Graves.] 

110 

SOLVE the triangle dC; given a, A and the product 7 of the bisectrices 

drawn from 2, C to CA, AB. 
SOLUTION. 
With the customary notations we get from the triangles BCE, BCF 


3 sin C r sin B 





a sin(d + LB) ‘ @- sin (4+ 4C)’ 





























SOLUTIONS OF EXERCISES 


pr sin B sin C 
a ~ sin(d + 48). sin(d + 4C): 


and therefore 


4%, _ cos? 4 (4 — C) — cos? 4 (B+ C) 
227 cos $(2 — C) + sin 3d 


Conversely from this relation } (4 — C) can be computed, and as 4 (4 -+ C) is 
also known, # and C, and from them the other parts of the triangle, are deter- 
mined. [Henry Heaton; T. U. Taylor.| 
112 
GENERALIZATION, 
To inscribe in the given triangle dC a triangle PGR whose sides OR, RP, 
PQ make given angles with BC, Cd, AL. 
SOLUTION I. 


The conditions of the problem give the directions of the sides of POR. Draw 
a parallel to AP, intersecting AB in P’ and CA in XR’. Through 7” and A’ draw 
parallels to OX and PQ, intersecting at VU’. Through Q’ draw a parallel to BC, 
intersecting CA at C’. Through C draw a parallel to C’?’. This parallel will in 
tersect AP at P. The remainder of the solution follows immediately. 
[A. B. Evans; E. L. Stabler ; Ormond Stoni a 


SOLUTION II. 


Draw the triangle ?’Q’R’ whose sides Q’R’, R’?’, P’Q’ pass through A, 2, ¢ 
and are parallel to OR, RP, PQ, and another 4’2’C"’ whose sides 2'C’, C’A’, 
A'B' pass through ?’, Q’, X’ and are parallel to BC, CA, AB. Divide the sides 
of ABC at P, Q, R in the same ratios in which the corresponding sides of A’A’C’ 


>=? 


are divided by ?’, Q’, RX’. PQR is the triangle sought. [Sarah Szold.] 
113 
THE transversal .1/V meets the sides A/, AC of the fixed triangle A?C, and 
makes 
AM .AN=BM.CN. 
Find its envelope. [O. Root, Jr.] 


SOLUTION I 


Ad W) 7 
= oy) , the point .J/ divides the 


segment AF in the same ratio in which WV divides the segment C4; in other 


According to the given condition 


‘ words, the variable points 17 and NV form similar projective ranges on the fixed 


lines AB, CA. Hence, the variable line A/V which connects corresponding 
points of these ranges envelopes a parabola. (See, for instance, L. Cremona’s 
Projective Geometry, translated by Leudesdorf, Oxford, 1885, p. 128.) The lines 
AB and AC are tangent to this parabola in / and C respectively, these being the 








eth ee 

















SOLUTIONS OF EXERCISES. 


g2 
points corresponding to the point of intersection A of the two ranges. 
[Alexander Ziwet.] 


SOLUTION II. 


Referring all parts of the figure to AZ, AC as co-ordinate axes, and putting 
AB=6, AC =<, 


AN= 1, 


AM= im, 


we have for the equation to 7V 
we 


and the parameters 7, 7 are connected by the relation, 


MW i 
= I. 


bt 
Differentiation and elimination give 
mm w- fix, oy. 
bx cy NO Ne’ 
whence for the equation to the enveloping parabola we have 
lar rs | V a 
Ne 7 VC FT [W. M. Thornton. 
114 
A CIRCLE intersects a conic in four points, ?,, 2, ?;, 2. Show that if the 
v-axis be parallel to the axis of the conic the area of their quadrilateral is 
(a, — x,) (1 —J’s)- [R. 77. Graves.] 
SOLUTION. 
Let the +-axis be parallel to, or-coincident with, the axis of the conic, and 
let the co-ordinates of the vertices taken in order be (44, 7), (4, 3’2), (4% 1g), (4 J’g)- 
Let 4 be the acute angle that each diagonal makes with the +-axis.* Then the 
area of the quadrilateral = } product of its diagonals into sin 24 
pa Se AO de ots 
— = 9 con® ae = (4% — 4) (1 — J). 
The same formula may be proved otherwise for the parabola as follows : — 
Then 


Let »* = 4ax be the equation to the parabola. 


Oo I O I O I Oo I 
I o 41 oOo lt #2 S-@ 
Area=+ 4 =H te 
, , a , , 
Ii V2 Vs Is Ky V2 Ix I 
. : : Ce ee he 
A 4p 43 ON J Ja SI 4 


*Salmon’s Conic Sections, 244.—[W. M/. 7:] 



































SOLUTIONS OF EXERCISES. 


oO oO o 1 

1| oO Oo 1 oO} } 1 I 
— oan 8a (i— = Sn (J — Js) (2 —J;,) 

. IIs V2 Is Is OSs | F ra +I3 Vo tI 

IIs I-Ie Is Ie | 

I | 2 I 
=+,| (W1 — Is) (J2 —I 

8a INtFIVet Ist, Wet ' ‘ 

I I 
= + —(y¥, — Is) (12 —%) 

- oO Ist Is 

I 
—=—+ = (ni — Is) (2? — 972) = & (4, — 4) 1, — J): [R. H. Graves.] 
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Suow that the angle between the tangent to the ellipse xr? + a*7* = a*/* 
and the tangent at the corresponding point of the principal circle is given by the 
formula, 

atan ¢ 
(1 —a) + tan? ¢’ 





tand = 


where 4 = (a — A) a is the ellipticity, and ¢-is the eccentric anomaly. 
SOLUTION, 


Let (2’, 9’) be the co-ordinates of the point of tangency with the principal 
circle; (2’, 7’’) those of the corresponding point of the ellipse. The equations 
of the tangents to the circle and the ellipse are 


7 1 a 
J oar yr" 1 Ty ‘ 
ifr! 6° 
and IA—ayrtt+s- 
ay y 


Hence the tangents /, g of the angles which the tangents to the circle and the 
ellipse make with the axis of 4, are 





a’ fx! 
P= y’ ’ = =< ay!" . 
But a’=acosg, jy’ =asing, 7’ =dbsing; 
b 
=— », and g=———. 
tan ¢ atang¢g 
g—fp atan¢ 





Then tand = 


I +pgo (1 —a) + tan? ¢” 
[William Hoover ; C.D. Schmitt. ] 

















94 SOLUTIONS OF EXERCISES. 


118 
SHow that the eccentric anomaly which gives the maximum deviation is 
tan~', (1 — 4), and find the co-ordinates of the point of contact 7 and the equa- 
tion to the tangent 7. 
SOLUTION. 
The deviation will be a maximum when tan ¢@ is a maximum. Equating the 
first differential coefficient to zero, 





a{(1 — a) + tan’ ¢] sec* ¢ ree .2tan¢ sec’ yg = 








whence ¢ =tan~'; (1—a@); 
n JE I — ‘|. I 
Ss i= ‘c= =, 
+ — th Vy (2 — a)’ 
. a i a(i— a) 
and v= ————_. , J na gt SIONS ie: 
y (2—4) Vv (2—4) 
The equation of the tangent ¢ is therefore 
a(i— a)? ale . 2h2 
, ————. += a’. 
) (2—a) JT ) (2 —a) ; [William Hoover. | 
119 
SHow that the portion of ¢ intercepted between the axes equals in length 
the sum of the semi-axes a, 4, and is divided at 7 into these two parts. ~ | 





SOLUTION. 


Let x =o and 7 = O in order in the equation of the tangent. The correspond- 
° l (2 tents “) 5 Tl / 2 9 
—,%=ay(2—4). ten = 7/ (4 + 0 
a ae SAV (2—M). Then t= (nF + 9) 
=a-+%, If <bethe segment adjacent to the axis of +, s: y’’ :: t: 3, whence 
> =~; and the other is a. [Wilham Hoover.| 
120 
SHow that the circumcircle of 4vy passes through the centre of curvature at 
7, and that the area of the circle of curvature equals that of the ellipse. 





ing intercepts are ), = 


SOLUTION, 
If » =the radius of curvature, we have for the ellipse, 


(a sin’ ¢ + 4 cos’ ¢)! 
= Pe == J (ad), 





by using sin ¢, cos ¢ given in solution of 118. The radius of the circumcircle 
of try = 4 a (a + 6), and the distance of its centre from 7=4(a@a— 6). .°. since 
plat 3), 4 4 (a@-— 6) are the sides of a right triangle, the circumcircle passes 
through the centre of curvature 7. We have zy? = zad. [William Hoover.] 





























EXERCISES. 


121 
SHow that when the ellipse varies subject to the condition, a + 4 constant, 
the locus of 7 is the hypocycloid, 


xri+ pi =(a + O)i, 
which touches each ellipse at the point 7. 
SOLUTION. 


2 y2 
The equation of the ellipse is * - 4 == I, (1) 


and we have a+6=c=aconstant. (2) 


Differentiating (1) and (2), a and 4 only variable, 


x ay” 

2 da=— B db, (3) 
da = — db; (4) 
2 apr ; 
hae ot (5) 


2 
(1) is 5 (a+6)=1, or b=ctyi, and a= chai. 


These in (2) give +4 + 7 =(a + +)', a hypocycloid. From this equation 
dy yi 


’ a a(i —a) 


= - thich for « = 2? = —— pes Of es epee Cc _ 
P a whic r )(2—4) a.) J (2 a) , becomes 
— , (1 — 4), showing that the ellipse and hypocycloid have a common tangent 
at 7. [Welham Hoover.] 


EXERCISES. 





137 
INTEGRATE the differential, 
_ # cos i 
sint # + cos‘ # [V. Root, Jr.] 
138 
Asout the vertices of an equilateral triangle three spheres are drawn with 


radii equal to the side of the triangle. Find the volume common to them all. 
[W. M. Thoruton.] 














96 EXERCISES. 


139 
AN ellipse touches the given ellipse a®y? + @1? — a’@ at the extremities of 
a diameter and has its focion the curve. Find the position of its axis major, and 
the position of the foci when the angle between the given diameter and the axis 
major of the given ellipse is a maximum. [R. H. Graves.] 
140 
In the triangle A/C, the radius of the inscribed circle is 
1(6+¢) cos’ 4A + (c + a) cos’? $8 + (a + 6) cos? $C 
3 cot hd + cot $42 + cot $C i 


[Ormond Stone.] 
141 


SOLVE the differential equation, 
ud*y + 2dudy=vdx dy, 


u,v, ¥ being all functions of +. [Wilham Hoover] 
142 
Let f be a prime number of the form 4z — 1, 7 being an integer. Let the 
primitive p” roots of unity be 





voy apy agyA? opr? — 
a a ea a i a 





Put # = 27) f, and let the negative terms of the series, 


sinf, sin(##), sin(/*#), 2. . sin(A?-*A), 
ew’. oe". a". ...3 Foe 
p=[tant(p+1)0—4g(e+oa"4+a0"+.. .)f 


Verify this relation for the particular cases f = 3, 7, 11, and give a general proof. 
[George Paxton Young.] 
145 
Ir a curve of the fourth order have a singular point at which are three coin- 
cident tangents, this line passes through the intersection point of the bitangent 
with the line of junction of the two points of inflexion of the curve. 


[¥. H. Loud.] 





144 


Suow that the diurnal path of the shadow of the top of a vertical rod ona 
horizontal plane will be a rectangular hyperbola, if 


2sind=sin|~+¢ 7; |. 
[E+¢)}+s (: | 


where ¢ is the latitude of the place and 0 the declination of the sun, assumed to 
be constant. [R. H. Graves.] 


bl ae ell 
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